Bond electronic polarization induced by spin 
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We study theoretically the electric polarization induced by noncollinear spin configurations in the 
limit of strong Hund coupling. We employ a model of two magnetic ions sandwitching an oxygen 
ion. It is shown that there appears a longitudinal polarization P x along the bond, which is roughly 
proportional to (m*) 2 — {mf) 2 where m*n) is the a;-component of the spin orientation vector at 
right (left) magnetic ion. A numerical study of the model Hamiltonian yields both longitudinal and 
transverse electric dipole moments. The transverse polarization is shown to have a non-uniform 
as well as a uniform component, with the latter being consistent with the previous theory. The 
longitudinal polarization is non-uniform and oscillating with the period half that of the spin order, 
but the local magnitude is typically much larger than the uniform transverse polarization and may 
be detected by X-ray/neutron scattering experiments. 

PACS numbers: 75.80.+q, 71.70.Ej, 77.80.-c 
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I. INTRODUCTION 

A number of recent experimental breakthroughs has 
revived interests in the phenomena of coupling of mag- 
netic and electric (dipolar) degrees of freedom in a class of 
materials known as "multiferroics" 1-7 . Some noteworthy 
observations include the development of dipole moments 
accompanying the collinear-to-helical spin ordering 3 ' 4,6,7 
and adiabatic control of dipole moments through sweep- 
ing of applied magnetic fields 1,2 , which all unambigu- 
ously point to the strong coupling of electric and mag- 
netic degrees of freedom in these compounds. A number 
of phenomcnological 8,9 and microscopic 10,11 theories has 
been advanced to establish the connection between non- 
collinear spin order and ferroelectricity. 

In particular the work of Katsura, Nagaosa, and Bal- 
atsky (KNB) 10 proposed a microscopic theory for the 
interplay between non-collinear magnetic order and the 
dipolar polarization of the electronic wave function in- 
duced by it. The magnetic (M) ion is modeled by three 
degenerate ti g levels experiencing some external mag- 
netic field (to guarantee magnetic order) and subject to 
spin-orbit coupling. Two such magnetic ions are bridged 
by an intermediate oxygen (O) atom which itself has no 
spin-orbit interaction. Solving the model Hamiltonian 
perturbatively in the M-0 hybridization amplitude, KNB 
finds an electronic polarization orthogonal to the M-O-M 
axis in the ground states of one and two holes. 

In this paper we re-visit the M-O-M cluster model of 
KNB, but in the different limit of a strong Hund cou- 
pling. The M-O-M model in this limit is exactly solvable 
in the absence of spin-orbit coupling, with two classes of 
degenerate eigenstates. For these states there is no net 



electric polarization. The spin-orbit coupling on the mag- 
netic sites is then introduced as a perturbation within 
each degenerate manifold. The problem is exactly solv- 
able again and we can calculate the polarization (r) for 
each of the eigenstates thus obtained. It is shown that a 
non-zero polarization develops along the direction of the 
M-O-M cluster. We denote such polarization as "lon- 
gitudinal", to distinguish it from the "transverse" po- 
larization, perpendicular to the cluster axis, obtained in 
previous theories 8,10,11 . 

Numerical study of the mean-field Hamiltonian for the 
cluster confirms the existence of longitudinal polarization 
predicted analytically, and reveals a non-uniform compo- 
nent in the transverse polarization which is unexpected in 
previous theories 8,10,11 . Namely, oscillating components 
appear in both the longitudinal and transverse polariza- 
tions with a vanishing macroscopic average whereas the 
polarization predicted by existing theories gives a "uni- 
form" component induced by a non-collinear spin order. 

This paper is organized as follows. In Sec. II, we 
introduce the mean-field Hamiltonian for the three-site 
cluster, and in the limit of a large Hund coupling, 
its approximate eigenstates are derived analytically in 
a perturbation theory with respect to the spin-orbit 
interaction. Electronic polarization is then calculated 
for these eigenstates to the lowest order in the spin-orbit 
coupling. In Sec. Ill, numerical solution of the mean- 
field Hamiltonian is presented. Fourier analysis of the 
polarization reveals the existence of both uniform and 
non-uniform components in the transverse part, and 
non-uniform component in the longitudinal part. We 
conclude with a discussion in Sec. IV. Some technical 
details of the derivations used in Sec. II are described in 
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FIG. 1: Arrangement of two magnetic ions and an oxygen 
atom in the M-O-M cluster, rhi and rh r are the magnetic 
moment orientations of the magnetic sites bridged to the in- 
tervening oxygen atom. Polarization P x along the cluster's 
axis is found (see text for details). 

the Appendix. 



II. MODEL AND ANALYTIC RESULTS 

As a minimal unit giving rise to spin-polarization cou- 
pling KNB introduced the model shown in Fig. 1. To 
guarantee that a net average magnetic moment exists, 
each magnetic site is subject to a spin polarization field 10 



H M = -UY J ™a-\ ]T S a , b (1) 
a—l,r \b—xy,yz,zx J 

where a = l,r refers to left and right magnetic site in 
the M-O-M cluster as depicted in Fig. 1. Spin-orbit 
coupling is assumed to exist for the magnetic ions which 
are connected to the oxygen through the wave function 
overlap. The spins S a ^ in each of the three degenerate 
t 2g orbitals are subject to the local magnetic field along 
in a = (sin0 a cos<j) a ,sin9 a sm(j) a ,cos8a)- Ultimately such 
effective fields arise due to inter-atomic exchange interac- 
tion but also due to the strong Hund coupling within the 
individual magnetic ion. In such a case the strength of U 
mimicking the Hund coupling energy may well exceed the 
spin-orbit interaction strength A. The large-JJ limit thus 
provides a natural starting point for the analysis of the 
spin-polarization coupling in the cluster model. In this 
section we propose a method that is particularly suited 
to treat the large U limit, namely U ^> A, V where V 
represents the M-0 hybridization. 

The overall Hamiltonian describing the three-atom 
cluster is given by H = Hm + Ho + Hy\ Hm is al- 
ready given in Eq. (1) and the other two terms, referring 
to the oxygen orbital and the M-0 hybridization, are 

b=x,y,z cr 



(7 

~ d^xyaPyv ~ < ^r,zxa'Pza + h.C.). 

(2) 

Particle-hole transformations have been implemented on 
all three atomic sites 10 . The on-site energy of the oxygen 
atom, — E p , becomes +E p after the transformation. The 
spin-orbit interaction, not included here, will be intro- 
duced later as a perturbation. 

Now we discuss how one obtains a low-energy effective 
Hamiltonian from Eq. (2). First, due to the large energy 
difference 2U separating the spin- up (parallel to rh a ) and 
spin-down (antiparallel to rh a ) t 2g states, we can truncate 
the high-energy d-orbital states from the outset and write 
down the effective Hamiltonian using only the low-energy 
operators 

dt.xy = COS y d+ xy] + e"**" sin y d+ xyL 

<» = cos y<^T + e ~ 40asin y d t^i- ( 3 ) 

The d yz orbital does not hybridize with any of the oxy- 
gen p-orbitals in the linear geometry we consider in this 
paper. This, and the p^-orbital which does not hybridize 
with any of the c?-orbitals, will be left out. 

In this reduced Hilbert space the low-energy effective 
Hamiltonian becomes 

H' = H' M + H'q + Hy 

— (da,xyda,xy + d^ zx d a ,zx) 

a—l,r 

H' = E p J2 {VyoVyv + PtaPza) 
a 

H 'v = V(d+ xyPl , y + d+ zxPhz 

~ <£,xyPr,v ~ d t,zxPr,z + h.C.). 

(4) 

In writing down H v spin rotations similar to Eq. (3) 
have been employed also for the p-orbital operators: 

4- ®a 4- —idi • 4- 

Pa,y = COS Y P yT +€ Sm y P 2/i 

PXz = cos—p^+e v-sm—p^. (5) 

The axis of the three-atom cluster is taken as the +x axis 
going from left to right in Fig. 1. Until the spin-orbit 
coupling is introduced, the (d xy ,p y ) orbital pair remains 
decoupled from (d zx ,p z ). 

The p-orbital operators {pi, y ,p r ,y) are not orthogonal 
to each other. The same is true of {pi, z ,p r , z )- For further 
manipulation of the Hamiltonian we need to introduce a 
set of orthogonal operators 
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(7) 

for d-orbital states. We have defined e 1TI — k/\k\ where 

K = (Pl,y\Pr,y) = (Pl,z\Pr,z) 

= cos^cos^+e^-«sin^sin^. (8) 
2 2 2 2 w 

Using these new operators, the Hamiltonian in Eq. (4) 
becomes 



H — Hy + Hz 

H y = - E (^U D y xv 



a=l,2 



H a 



E ( ^ z « ) w « 

Q = l,2 



Z n 



= Eq + E a 



cos /3 Q — sin fi a 
— sin Pa — cos /3 a 



(9) 



with E = (U - E p )/2, E a = y/(U + E p y 2 /A + V*, 
and Vi, 2 = Vy/l±\K\, (cos/? Q ,sin/3 Q ) = {(U + 
Ep)/2Ea,Va/E a )- Hy and Hz can be diagonalizcd 
through the rotation 



D 



ct.xy 



D, 



a,zx 

Zry 



— cos/3 Q /2 sin/3 Q /2 \ 
sin/3 Q /2 cos/3 Q /2 y 

— cos /3 Q /2 sin/3 Q /2\ 
sin/3 Q /2 cos/3 Q /2 y 



0a 

Va 
Va 



(10) 



=]T-(£o+£a)^a+E(^- £; °)^«' 
a a 

ffz = E -(Eo + EcWWa + E(^- 

(ii) 

Four degenerate levels (-01 , i>2, ip[ , ip'2) are separated from 
the other four degenerate set (ipi,tp2, fi, ¥2) w ^h an en- 
ergy spacing of nearly U + E p . This large energy separa- 
tion sets the stage for introducing spin-orbit interaction 
Hso = XS • L within each of the four-dimensional mani- 
folds, but not between the two manifolds. We obtain the 
effective Hamiltonian valid within each manifold, 



H 



02 + 
V>i + 
\^2 + J 

+ 
V2 

\^2 ) 



( -Ui -iAi -i\ 2 

— U2 — iX 2 — iA 3 

i\i iX 2 —U\ 

V i\ 2 iA 3 -U 2 



02 
V>i 
02 



(E\ 


iAi »*2 s i 

V a' ix' 



i — iA 2 





El -iX' 2 -iX' 3 



E P 2 



<Pi 
¥2 

¥2 



(12) 



where U a = E a + E a , E v a = E a - E , and 



iX 1 = M |if S o|0i) = y cos 2 ^( m x + m?) 

iX 2 = {^Hsofo) = mHsoWi) 
iX fix fi 2 , x x , 

= ~2 cos T cos T ^ mr ~ mi ' 

W 2 \H SO \rk) = y cos 2 ^ (m* + mf) 



iA 3 = 

iAi = 
iA' 2 = 



iA', = 



^i|i/ 50 |^) = ysin 2 ^(m-+mf) 
(yi I -ffso |<^2) = (<p2\Hso\<Pi) 

mf) 



iX . fix , fi 2 , a 
— sin — sin — (m r 



(^ 2 |iJ S o|<^2) = y sin 2 -j (m x r + mf 



(13) 



Here m x a refers to the x-component of the local quanti- 
zation axis for left and right magnetic sites. Because of 
small V/U, we can safely take 



Ex = E 2 =E=(U + E p )/2 

COS fia/2 W 1 

sinfia/2 w ^/(IZ + Bp) 



(14) 



in Eq. (13). The above Hamiltonian H can be diagonal- 
izcd to the form 



into 



h = y j ^- + E E " ■ 



(15) 
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with the eigen-states 



*i - ^1+^2 + ^1+^2 

One can verify that the coefficients satisfy 



(16) 



Kii = ia if 2 i = isi. 
\a\ = \a\, idi = idi. 



Cid = CkleRe. (17) 



This completes the formal derivation of the full set of 
eigenstates in the limit U 3> A, V. The remaining task 
is to calculate the polarization (r) for each of the eigen- 
states obtained. 

Due to the shapes of d- and p-orbitals, only the follow- 
ing overlap integrals are non-zero (x = xx): 



{dr : xya \^\Pya) — (dl : xya \X-\Pycr) 
= (d r ,zxa\x\Pza) = {dl, zx<T \x\p z<T ) = Lx. (18) 

Here L = J d 3 rd ijXya (r)xLpy a (r) (i = r,l). If non-zero po- 
larization develops, it can only be in the x direction, i.e. 
along the axis of the three-atom cluster, within the low- 
est order in A/A, A = U + E p . Note that the transverse 
component obtained in Ref. 10 corresponds to the terms 
of the first order in A/ A. This contribution is numeri- 
cally studied later in Sec. III. Analytic calculations of the 
diagonal and off-diagonal contributions to the polariza- 
tion are given in Appendices. The final expression for the 
polarization P reads (details are given in the Appendix) 



P* ( = (*i|r|* i ) = 2^[<^i|r|^) + ^i|r|^>] 



*2 = (»V'l-*^-'0i+^2)/2, E^ = 

*3 = (-*^i + #2 - $ + V> 2 )/2, E\ = 
* 4 = HV'i - #2 + + V> 2 )/2, El = 



-U-X 1 -X 2 

-u+x 1 -x 2 
-C/-A1+A2 
-C/+A1+A2. 

(21) 



For each of the eigenstates above the polarization Py i 
reads 



P*X = P* 4 



2LU 



a; x 



P\I>2 — P 1 !'^ — X X 



A 

2lf 



(22) 



The polarization, constant in magnitude but reversible 
in sign, is developed along the M-O-M cluster. Since Ai 
and A 2 depend on the spin directions and in particular 
on mf and mf as shown in Eq. (20), which one of the 
four states shown in Eq. (21) becomes the true ground 
state depends on the set of values (mf, mf). In Fig. 2 we 
map out the polarization directions of the ground state 
for all the available situations within the zeroth order 
in A/ A. The reversal of the polarization direction occurs 
whenever \mf\ = \mf\ line is crossed. The sudden change 
of the sign is due to level crossing of the two eigenstates 
when Ai or A2 vanishes at \mf\ — \mf\. The ground 
state polarization is consistent with the functional form 
P x <~ (mff - (mf f. 

If a single hole is introduced in the cluster, the po- 
larization is given by one of the expressions in Eq. (22). 
With two holes, two opposite polarizations from Eq. (22) 
cancel out to produce P = 0. In the numerical results 
presented in the next section, however, one finds a fi- 
nite longitudinal polarization even for the two-hole case 
because of the higher-order terms in A/A, and the polar- 
ization direction is still dictated by P x ~ (mf) 2 — (mf) 2 . 



„ . fx ^£^5. (i9) 



Note that and C1C2 are a l so dependent on the spin 
configuration. The full spin dependence of the polariza- 
tion may be rather complicated. However, because of 
small V/U, we can take U\ = U 2 = U and 
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l\ 1 




-1 I /' 
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/ + 

k. 


+ \ 






-1 



mf 



Ai 



^3 = -(mf + mf) 



A 2 = -(mf-mf). 



(20) 



Under such an approximation the low-energy eigenstates 
are given by the surprisingly simple form 



FIG. 2: The direction of longitudinal polarization P x = P • x 
in the ground state wave function plotted for (mf,mf), — 1 < 
mf,mf < 1. +(-) refers to P x > (<) 0. P x is indeterminate 
along two dotted lines (mf = mf and mf — —mf) due to the 
level crossing occurring along these lines. 
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III. NUMERICAL RESULTS 



The results of the previous section gave the electronic 
polarization solely along the axis of the atomic cluster. 
In the analysis it was assumed that the large effective 
Zeeman energy U dominates over both the spin-orbit in- 
teraction A and the hybridization amplitude V. It was 
also assumed that d yz and p x orbitals did not play a 
role. According to Ref. 10, on the other hand, trans- 
verse polarization such as P y can only arise from the 
overlap integral (d yz \y\p z ) or (d xy \y\p x ). In the numer- 
ical approach this condition can be relaxed and indeed 
we find non-zero transverse polarization as well as the 
longitudinal polarization in accordance with the theory 
presented in the previous section. In general the results 
depend sensitively on the number of holes inserted in the 
cluster, with the even and odd number of holes yielding 
different behavior. As a representative case for each, we 
consider one- and two-hole situations in detail. When ex- 
trapolated to the infinite lattice they correspond to the 
metallic and insulating cases, respectively. 

The Hilbert space for the single-particle mean-field 
Hamiltonian is 18-dimensional, consisting of Eq. (1) and 
Eq. (2), with six from each atomic site. In practice, one 
can treat the 16-dimensional problem leaving out the oxy- 
gen p x orbital without loss of generality, since it remains 
decoupled from the rest of the atomic states in the exact 
solution. 



A. One hole (double-exchange interaction) 



The polarization (r) is calculated with respect to the 
lowest-energy eigenstate of the particle-hole transformed 
Hamiltonian corresponding to putting one hole in the 
cluster. For ease of presentation we restrict ourselves to 
situations where both spins lie in the XY plane: rhi — 
(cos <f>i, sin^, 0), rh r — (cos </> r , sin0 r , 0). The resulting 
longitudinal polarization P x is plotted as a function of 
(<pi, 4> r ) in Fig. 3. The similarity between the numerically 
obtained P x and the predicted behavior P x ~ cos 2cp r — 
cos 20; is striking. The magnitude of the polarization, 
according to Eq. (22) and using the parameters used in 
the numerical study, is expected to be \P X \ ~ 0.05L, in 
good agreement with the numerical results. 

In addition, we find non-zero P y in the numerical calcu- 
lation due to the inclusion of d yz orbital. The behavior 
of P y we find is consistent with the directional depen- 
dence P y ~ sgn[cos(A0/2)] sin(A0/2) (A<f> = <f> r - </>;) 
expected in Ref. 10 (See Eq. (4) in Rcf. 10). Also the 
order of magnitude \P y \ ~ 0.01X obtained numerically is 
consistent with Ref. 10. P z vanishes for the two spins 
lying in the XY plane, as shown both analytically and 
numerically. 




FIG. 3: (color online) The calculated polarization P x and P y 
from exact diagonalization of the full M-O-M cluster Hamil- 
tonian with a single hole. Dotted white lines indicate P x 
or P y equal to zero. P x vanishes not only when the spins 
are collinear, but also when <f>i + <j) r = ±7r, as discussed in 
Sec. II. The overall directional dependence for P x is consis- 
tent with cos 2(j> r — cos 2(f>i . The sharp change in the sign 
of P x along 4> r — <pi ± 7r lines is due to the energy level 
crossing in the ground state. P x turns to zero along the 
lines (indicated as dotted white lines) where the spins be- 
come collinear. The parameters used in this calculation are 
(U/V, E p /V, X/V) = (10, 30, 0.3). 



B. Two holes (superexchange interaction) 

When two holes are introduced in the cluster, the 
ground state can be changed adiabatically without a 
level crossing as the two spin directions are varied, and 
hence the polarization exhibits a smooth change. Typi- 
cal numerical results for the longitudinal and transverse 
components of the polarization are shown in the upper 
and lower panels of Fig. 4. As in the previous con- 
sideration for one hole, the spin directions are varied 
within the XY plane with rhi = (cos <pi , sin <f>i , 0) and 
fh r = (cos 4> r , sin 4> r , 0) . 

We analyzed the directional dependence of the polar- 
ization by using the following forms: 



P x /L = v4[cos(20 r ) -cos(20 z )] 
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FIG. 4: (color online) Induced polarization P x and P y in a 
unit of L for the superexchange case with two holes in the 
three-atom cluster. The parameters used in this calculation 
are the same as in Fig. 3. 

Py/L = -B 1 sin((j> r -<f>i)+B2[sm(2(f> r )-sm(2<f>i)]. 

(23) 

While the least-square fit to the numerical results were 
made using as many as six different basis functions, only 
those listed above made a sizable contribution. The angle 
dependence of the transverse polarization P y predicted in 
Rcf. 10 corresponds to the B\ sin(</> r ~ 4>i) term above, 
whereas the angle dependence of P x is consistent with 
the single-hole results analytically derived in the previ- 
ous section. The second term in P y , proportional to B 2 , 
was not anticipated in analytic theories but show up in 
the numerical analysis. In the spiral phase for the spins, 
Bi sin(0,. — 4>{) produces a uniform transverse component 
while all other terms give rise to a macroscopically van- 
ishing polarization. In this regard the new terms (A and 
B 2 ) correspond to "oscillating" component of the polar- 
ization in contrast to the uniform component given by 
Bl 

Shown in Fig. 5 are the coefficients, A, Bi, and 
B 2 , for a wide range of parameters, < X/V < 10, 
< U/V < 10, obtained by fitting the numerical data 
using Eq. (23). Several features are worth mentioning: 
(i) A and B 2 are remarkably similar over most of the 
parameter space. This implies that the non-uniform po- 
larization develops with similar magnitudes along and 
orthogonal to the bond axis, (ii) Large non-uniform po- 



larization components are attained for large U/V and 
small X/V as indicated by the blue and violet regions in 
Fig. 5. On the other hand, large uniform polarization 
(Si) is achieved for large X/V and small U/V in agree- 
ment with Ref. 10, where A > U was implicitly assumed. 
To be concrete, the largest values for A, B\ and B 2 are 
obtained for [U/V, X/V) = (1.795,0.256), (0.256,0.256), 
and (1.794,0.256) respectively, with the magnitudes in- 
dicated on the vertical bar next to each figure. Using 
an earlier estimate 10 , these maximal values correspond 
to A ~ B 2 ~ 2 x 10 2 nC/cm 2 and B x ~ 80 nC/cm 2 . For 
U /V large, the polarization mainly possesses the oscillat- 
ing components with a much smaller uniform component 
as explicitly demonstrated at the bottom of Fig. 5 for 
U /V fixed at a large value (U/V — 10) and varying X/V. 

The macroscopic polarization is governed by the B\ 
term alone. On the other hand, the Fourier analysis re- 
veals the local magnitudes of the oscillating component to 
be much larger than the uniform polarization over most 
of the parameter ranges, increasing the chance to ob- 
serve coupling of the local dipole order to external probes 
such as X-ray or neutron scattering. Particularly in spi- 
ral magnets with ordering wavevector Q, the oscillation 
of electric dipole moments appears with the wavevector 
2Q, since we have found that the nonuniform part has the 
form sin 2<f>i — sin 2<fi r for the transverse component and 
cos 2<fii — cos 2<fi r for the longitudinal one. Figure 6 shows 
the simplest example with the periodicity of the eight M 
ions where relative angles between the neighboring spins 
are 7r/4. In this case, electric dipole moments induced 
at oxyen ions also helically rotate with half the magnetic 
periodicity. Note that the amplitude of this nonuniform 
part is much larger than the uniform part. It is expected 
that oxygen ions subject to such electric dipole field shit 
in either parallel or anti-parallel directions to the local 
dipole moment. 



IV. DISCUSSION 

In this paper we have developed a theory for the spin- 
induced electric polarization in a simple three-atom clus- 
ter model proposed in Ref. 10. The truncation to the 
low-energy space is allowed by the large energy separa- 
tion between spin-up and spin-down states of the mag- 
netic site due to the strong effective Zeeman field, and 
the much smaller energy scale for the pd-hybridization 
V and spin-orbit interaction A. The net electric dipole 
moment is found to be along the cluster's axis within our 
approximation, and the numerical calculation confirmed 
this conclusion. The longitudinal polarization has an os- 
cillating character, and vanishes in the macroscopic limit. 
Numerical analysis further identified uniform and non- 
uniform components in the transverse polarization. The 
magnitude of the non-uniform transverse polarization is 
comparable to the non-uniform longitudinal one. Such 
non-uniform polarizations generically have larger magni- 
tudes on the atomic scale than the uniform transverse 
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component. 

In this paper we have concentrated on the t 2g electron 
systems because the spin-orbit interaction has the ma- 
trix elements within the ^-orbitals, while they are zero 
among the e s -orbitals. The effective Zeeman field at each 
site corresponds to the mean field from the surrounding 
spins, which leads to the similar model to the double ex- 
change model discussed in this paper. Generalizations to 
the more realistic models including the hybridization be- 
tween e g - and i 2g -orbitals and spin-orbit interaction on 
oxygen orbitals are now in progress, which will be applied 
to the realistic materials such as i?MnC>3 (i?=Tb,Gd). 
However, our present model will have relevance to the 
materials such as YVO3 12 , i? 2 Ti 2 7 (i?=Gd,Tb) 13 which 
have partially filled t 2g electrons. More serious consider- 
ations on each material including the finite bond angle 
etc. are left for future investigations. 

Recent experimental works on multiferroics with spi- 
ral spin structure with the ordering wavevector Q have 
reported the lattice modulation at the wavevector 2Q in 
addition to the macroscopic polarization 1 . Since dipolc 
moments are linearly coupled to the ionic displacements, 
it may correspond to the non-uniform dipole moments 
shown in this paper. In our scenario of magnetically in- 
duced dipole moments, the modulation of the oxygen ions 
is non-collinear and has also a spiral structure, reminis- 
cent of the MnSi case 14 . A similar lattice modulation but 
with the collinear shifts has been discussed by Sergienko 
and Dagotto 11 . Further X-ray and neutron scattering ex- 
periments are required to specify the detailed character 
of this lattice modulation. 

On the theory side, a Ginzburg-Landau consideration 
led to the conclusion that the uniform polarization direc- 
tion is orthogonal to the modulation wavevector of the 
spin 8 ' 9 . In Ref. 8 some non-uniform components of the 
polarization is predicted on the basis of Ginzburg-Landau 
theory without a detailed discussion of its implications. 
Furthermore, consideration of the orbital structure for 
the magnetic site is missing in the phenomenological the- 
ories advanced so far. For instance, the t 2ff levels alone 
break the rotational symmetry in the orbital space and 
the corresponding Ginzburg-Landau theory need not pos- 
sess the invariance under spatial rotation. It will be in- 
teresting to develop a Ginzburg-Landau theory with a 
proper symmetry consideration for the t 2g systems to see 
if the new features we observe, i.e. non-uniform longitu- 
dinal and transverse polarizations, emerge naturally. 

Our results may also open up a possibility for a 
nanoscale device of switching the polarization by means 
of an applied magnetic field or the magnetization by 
means of an applied electric field. The so-called "quan- 
tum magnet" might be a candidate to such systems. The 
present results even pave a way to fabricate crystaline 
magnets, ferroelectrics and multiferroics having the lon- 
gitudinal polarization. 
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APPENDIX: CALCULATION OF THE DIPOLE 
MOMENTS 

In this appendix, we give details of analytic calcula- 
tions of dipole moments associated with the approximate 
eigenstates of the total Hamiltonian. 



1. Diagonal terms 

Since 

(YiNA,*,) = - -it— [1 - 1 + \k\ - \k\\ x = 0, 
2^1 + |k| 

(A.l) 

we have 

(ViMVi) = [(Di^m) + (Fi|r|£> liX2/ )] = 0, 

(piM^i) = KDi^Myi) + (YilriLVj,)] = o. 

(A.2) 

All other diagonal terms vanish as well, (ipm\i"\ipm) — 

Wm\ r Wm) = {<Pm\r-\Vm) = (^4| r l^m) = 0. 



2. Off-diagonal terms 

In view of the symmetry of the orbital overlaps 
given in Eq. (18), the following averages must be zero: 
(Vi|r|Vi> = (V-i|r|V 2 > = (WWi) = (^|r|^ 2 ) = 0. 



3. Dipole moments for the eigenstates 

Thus, only (ipi\r\tp 2 ), (^il r IV> 2 ) and ((pi\r\ip 2 ), 
(<Pi\r\<p' 2 } may contribute to the polarization. Since 



(D hxy \r\Y 2 ) =-xx Ly/1-\ K \ 
(Yi\t\D 2 , xv ) =-xx Ly/1 + \ K \, (A.3) 



we have 



(Zi|r|£> 2 ,«x) =-xx Ly/l + \k\ (A.5) 



<Vi|r|^2> = S ives 

01 . 02 /r . I i v \ ■ & ^/vl in \ 

- cos — sin y (Di )X j,|r|y 2 ) - sin y cos y |r|D 2 ,x3/) 

- * x — [(1 - \k\) + (1 + N)] =« — , <^i|r|-02> = 

01 ■ 02 /r * \ \ry \ -01 02 I rj I I n \ 

, M , -cos — sin— (D lzx r Z 2 ) -sin — cos— {Z x r £> 2 ,*z) 

(<Pi|r|<^ 2 ) = 2 2 2 2 



■ 01 02 it, I I -*r \ i A . & /v i m \ ~ T 

smy cosy(Di ia;y |r|y2) + cosy smy(Yi|r|D 2 ,xy) ~ 1 x 

. 2LV" /- . V2LV r ,,,,,, 



2Ly 



(A. 4) siny cosy(£>i )Za; |r|Z2) + cos y sin y (^i|r|£» 2j2X ) 



Similarly, „ V2LV 



x x — - — VI - o"i • 02. (A.6) 



<£>i, ZI |r|Z 2 ) = -i x Lv/1- |k| 
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FIG. 5: (color online) (top three figures) Fourier components 
A, B\, and B2, in the polarization P x /L and P y /L using the 
fitting formula, Eq. (23). (bottom) Fourier components with 
U/V = 10 and varying X/V. 
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FIG. 6: Spiral spin structure (red arrows) with an angle shift 
of 7r/4 between neighboring magnetic (M) ions and the in- 
duced dipole moments (parallel or anti-parallel to black ar- 
rows) at oxygen (O) sites with half the periodicity. Uniform 
component is not shown since it is much smaller than the 
nonuniform part. 



